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ABSTRACT 

In the context of the study of the properties of the mutual mass distribution of 
the bright and dark matter in elliptical galaxies, I present a family of two-component, 
spherical, self-consistent galaxy models, where one density distribution follows a 71 
profile, and the other a 72 profile [hereafter (71,72) models], with different total 
masses and "core" radii. A variable amount of (radial) orbital anisotropy is allowed 
in both components, following the Osipkov-Merritt parameterization. For these 
models, I derive analytically the necessary and sufficient conditions that the model 
parameters must satisfy in order to correspond to a physical system (the so-called 
model consistency). Moreover, the possibility of adding a black hole at the center of 
radially anisotropic 7 models is discussed, determining analytically a lower limit of the 
anisotropy radius as a function of 7. The analytical phase-space distribution function 
for (1,0) models is presented, together with the solution of the Jeans equations and 
the quantities entering the scalar virial theorem. It is proved that a globally isotropic 
7 = 1 component is consistent for any mass and core radius of the superimposed 7 = 
model; on the contrary, only a maximum value of the core radius is allowed for the 
7 = model when a 7 = 1 density distribution is added. The combined effects of mass 
concentration and orbital anisotropy are investigated, and an interesting behavior of 
the distribution function of the anisotropic 7 = component is found: there exists 
a region in the parameter space where a sufficient amount of anisotropy results in 
a consistent model, while the structurally identical but isotropic model would be 
inconsistent. 



Subject headings: galaxies: elliptical - stellar dynamics - dark matter - black holes 



1. Introduction 



When studying a stellar dynamical model (single or multi-component), the fact that the 
Jeans equations have a physically acceptable solution is not a sufficient criterion for the validity of 
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the model: the minimal requirement to be met by a physically acceptable model is the positivity 
of the phase-space distribution function (DF) of each distinct component. A model satisfying this 
essential requirement (which is much weaker than the stability of the model) is called a consistent 
model; moreover, when the total potential is determined by the model density components through 
the Poisson equation, the model is called self-consistent. In other words, a self-consistent model 
represents a physically acceptable, self-gravitating system. 

Two general strategies can be used to construct a (self) consistent model, or check whether 
a proposed model is (self) consistent: they are commonly referred to as the "/-to— p" and the 
"p— to— /" approaches ( [Binney and Tremaine 1987| Chap. 4, hereafter BT87). An example of the 
first approach is the extensive survey of self-consistent two-component spherical galaxy models 
carried out by Bertin and co-workers (Bertin, Saglia, and Stiavelli 1992), where they assume for 
the stellar and dark matter components two distribution functions of the foo form (and so positive 
by choice; [Bertin and Stiavelli 1984| ). The main problem with this approach is that generally the 
spatial density is not expressible in terms of simple or at least well known functions, and so only 
numerical investigations are usually feasible. 

In the second approach, the density distribution is given, and assumptions about the model 
internal dynamics are made, making the comparison with the data simpler. The difficulties 
inherent in the operation of recovering the DF in many cases prevent a simple consistency 
analysis. In particular, in order to recover the DF of spherical models with anisotropy, the 
Osipkov-Merritt technique ( Osipkov 1979] , Merritt 1985a , hereafter OM) has been developed 



from the original Eddington (1916) method for isotropic systems, and widely used. Examples of 
numerical application of the OM inversion to one and two-component spherical galaxies can be 
found in the literature (see, e.g., Ciotti and Pellegrini 1992, hereafter CP92; [Hiotelis 1994 ; Carollo, 



de Zeeuw, and van der Marel 1995a; Ciotti and Lanzoni 1997, hereafter CL97). If one is just 



interested in the (self) consistency of a stellar system the previous methods obviously give "too 
much", i.e., give the full DF. In the OM framework, a simpler approach in order to check the (self) 
consistency of spherically symmetric, multi-component models, is given by the method described 
in CP92. This method uses directly the radial density profile of each component and of the total 
potential, and gives necessary and sufficient conditions for the model (self) consistency, avoiding 
the necessity of recovering the DF itself. Moreover, since it requires only spatial differentiation 
and inequality checks, this method is best suited for analytical investigations. 

The importance of studying multi-component galaxy models cannot be underestimated: in 
fact it is now accepted that a fraction of the mass in galaxies and clusters of galaxies is made of a 
dark component, whose density distribution - albeit not well constrained by observations - differs 
from that of the visible one (see, e.g., Bertin et al. 1994, Carollo et al. 1995b| , [Buote and Canizares 



1997] [Gerhard et al. 1998| ). Moreover, there is increasing evidence for the presence of massive black 



holes (BHs) at the center of most (if not all) elliptical galaxies (see, e.g., Harms et al. 1994, 



van 



der Marel et al. 1997a, van der Marel, de Zeeuw, and Rix 1997b, |Richstone 1998 ). It follows that 



the obvious generalization of the one-component spherical models (the dynamicists zero-th order 
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approximation of real galaxies) is not only in the direction of the actively developed modeling of 
axisymmetric and triaxial systems (see, e.g., de Zeeuw (1996) for a review) but also in the direction 
of the construction of two-component analytical models, and in the study of their phase-space 
properties, a field far less developed. Prom this point of view the 1 st order approximation of real 
galaxies is the construction of analytical, spherically symmetric, (self) consistent two-component 
galaxy models. 

Unfortunately, few examples of two-component systems in which both the spatial density and 
the DF are analytically known are available, namely the very remarkable axisymmetric Binney- 
Evans model ( Binney 1981 , Evans 1993| ), and the spherically symmetric two-component Hernquist 



model (HH models, Ciotti 1996, hereafter C96), and so it would be particularly interesting to 
find other members of this exclusive club. In C96 the (successful) choice of the Hernquist density 
distribution ( Hernquist 1990| , hereafter H90) as building block for a two-component model with an 



analytical DF, was suggested by the extremely simple (and algebraic) expression of its potential 
as a function of radius. Moreover, the application of the CP92 method to HH models revealed 
itself as both simple and fruitful, explaining many properties of the derived DF. In this line, a 
natural and promising extension of the HH models is obtained by considering the wider family of 
spherically symmetric, two-component, anisotropic (71,72) models. This family of models is made 
by the superposition of two 7 models [see equation (7)] with different total masses, scale-lengths, 
and slopes 7. The mass concentration and amount of mass of the two distributions are described 
by four free parameters, and orbital anisotropy is allowed in both components, following the OM 
prescription. Because the Hernquist profile is obtained from 7 = 1, it will referred as a 7 = 1 
model, and so with the adopted nomenclature the HH models discussed in C96 will be called (1,1) 
models. Note that an increasing interest in (71,72) models is seen in simulation and observational 
works: e.g., Pellegrini and Ciotti (1998) used (2, 1) models in their numerical simulations of hot 
gas flows in X-ray emitting elliptical galaxies, Loewenstein and White (1999) used galaxy models 
similar to (1,1) models in the inner regions in order to observationally constraint the properties of 
dark matter halos around ellipticals. 

As expected, it is not possible to find the analytical DF for (71,72) models in the general 
case, but the fact that the DF of (1,0) models with OM anisotropy is completely expressible in 
analytical way, as proved here, is still of great interest. The study of (1,0) models is also useful for 
many other reasons: to provide an analytical DF for a two-component system for which the virial 
quantities and the analytical solution of the Jeans equation can be found explicitly; to arrange 
initial conditions for numerical simulations of two-component systems; to investigate the role of 
anisotropy and mass distribution of each component in determining the positivity of their DF. 
The availability of the analytical DF for a two-component stellar system allows us to arrange with 
great accuracy the initial conditions for numerical simulations aimed at investigating the stability 
of galaxy models in the presence of dark matter halos, or with a central BH. A work on stability 
of (1,0) models is in progress (Londrillo and Ciotti, in preparation). 



Strictly related to the last point above, is the trend shown in the numerical investigations 
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of two-component models described in CP92, i.e., the difficulty of consistently superimposing a 
centrally peaked distribution to a centrally flat one. CP92 showed numerically that King (1972) 
or quasi-isothermal density profiles can not be coupled to a de Vaucouleurs (1948) model, because 
their DFs run into negative values near the model center. On the contrary, the DF of the de 
Vaucouleurs component is qualitatively unaffected by the presence of centrally flat halos. From 
this point of view, the C96 work on (1,1) models is complementary to the investigation of CP92: 
in the (1,1) models the two density components are both centrally peaked, and their DF is positive 
(in the isotropic case) for all the possible choices of halo and galaxy masses and concentrations. 

The implications of these findings have been not explored sufficiently One could speculate 
that in the presence of a centrally peaked dark matter halo, King-like elliptical galaxies should be 
relatively rare, or that a galaxy with a central power-law density profile cannot have a dark halo 
too flat in the center. In fact, observations of the central surface brightness profiles of elliptical 



galaxies (see, e.g., Ferrarese et al. 1994, Jaffe et al. 1994; Moller, Stiavelli, M., and Zeilinger 1995, 



Lauer et al. 1995, Kormendy et al. 1995, Byun et al. 1996| ), and bulges of spirals (Carollo and 



|Stiavelli 1998 ), as well as high-resolution numerical simulations of dark matter halo formation 



(see, e.g., Dubinski and Carlberg 1991, White 1996; Navarro, Frenk, and White 1997) seem to 



point in this direction. 

In this paper, I further explore the trend that emerged in CP92 and in C96, determining the 
limits imposed by phase-space constraints (i.e., the DF positivity) on the parameters describing 
the (71,72) models and 7 models with a central BH [hereafter (7,BH) models]. I focus on the (1,0) 
models, in which one component is centrally peaked (7 = 1 density profile), and the other has a 
more flat core (7 = density profile). With the aid of the derived, analytical DFs, more stringent 
conclusions are then reached. 

In Section 2, I briefly review the technique developed in CP92 and applied in C96 to (1,1) 
models, and I formulate it in a way suitable for its application to the present problem. In Section 
3, the (71,72) models are introduced, as well as the CP92 method used to discuss the limits 
imposed on their parameters by requiring the positivity of the DF of the two components. In 
Section 4, the DF for the two components of (1,0) models are derived explicitly, and in Section 
5, the exact boundaries of the region of consistency in the parameter space are obtained and 
compared to those given in Section 3. Finally, in Section 6, the main results are summarized. In 
the Appendix, the analytical expressions for the velocity dispersion profiles of both the (1,0) model 
components and the virial quantities useful in applications, are derived in the general OM case. 



2. The Consistency of Multi Component Systems 

As outlined in the Introduction, a stellar system described as a sum of different density 
components is called consistent if each is non-negative over all the accessible phase-space; 
a consistent, self-gravitating system is called self-consistent. The technique developed in CP92 
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permits us to check whether the DF of a multi-component spherical system, where the orbital 
anisotropy of each component is modeled according to the OM parameterization, is positive, 
without calculating it effectively. In the OM formulation, the radially anisotropic case is obtained 
as a consequence of assuming / = f(Q) with: 



Q = £ 



2rV 



(1) 



where £ and L are respectively the relative energy and the angular momentum modulus per unit 
mass, r a is the so-called anisotropy radius, and f(Q) = for Q < 0. With this assumption, the 
models are characterized by radial anisotropy increasing with galactic radius, and in the limit 
as r a — ► oo, the velocity dispersion tensor becomes globally isotropic. For a multi-component 
spherical system, the simple relation between energy and angular momentum prescribed by 
equation (1) allows us to express the DF of the k-th component as: 



/k(Qk) 



l 



d 



dg k 



1 [Q* d 2 g k d^r 



V8tt 2 dQ k Jo d^ T VQk ~ *t V8ir 2 



5r T ' 



where 



Qk(r) 



1 + 



Pk(r), 



(2) 



(3) 



ak , 



* T (r) = Efc^k(^) is the total relative potential, Q k = £ - L 2 /2rl k , and < Q k < * T (0). The 
second equivalence in equation (2) holds for untruncated systems with finite total mass, as the 
models discussed here (see, e.g., BT87, p. 240). In C96, the original CP92 technique was applied 
to (1,1) models using the relative potential $> k of the investigated component as the independent 
variable, but here the radius r is found to be a more convenient variable. So, we have now 



Theorem: A necessary condition (NC) for the non-negativity of f k is: 

dg k (r) 



dr 



< 0, < r < co. 



If the NC is satisfied, a strong sufficient condition (SSC) for the non-negativity of f k is: 

dg k (r) rV* T (r) 



d_ 

dr 



> 0, < r < oo. 



dr Mx(r) 

Finally, a weak sufficient condition (WSCQ) for the non negativity of f k is: 



(4) 



(5) 



d_ 

dr 



dg k {r) r 2 
dr M T (r) 



> 0, < r < oo. 



(6) 



lr The WSC is better suited than the SSC for analytical investigations, due to the absence of the weighting square 
root of the total potential. 
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Proof: See CP92 and C96. 

Some considerations follow from looking at the above conditions. The first is that the violation 
of the NC [equation (4)] is connected only to the radial behavior of pk and the value of r a k, and 
so this condition applies independently of any other component added to the model. Obviously, 
the condition imposed by equation (4) is only necessary, so /k can be negative (and so the k 
component will be inconsistent) even for values of model parameters allowed by the NC. This 
is due to the radial behavior of the integrand in equation (2), which not only depends on the 
particular and r a k, but also on on the total potential. To summarize: a model failing the NC 
is certainly inconsistent, a model satisfying the NC can be consistent. The second consideration 
is that a model satisfying the WSC (or the more restrictive SSC) is certainly consistent, a 
model failing the WSC (SSC) can be consistent, due to the sufficiency of the conditions given by 
equations (5)- (6). As a consequence, the consistency of a model satisfying the NC and failing 
the WSC (or the SSC) can be proved only by direct inspection of its DF. For example, if one 
find that for r a k < ^(NC) the model is inconsistent [i.e., equation (4) is not verified], while 
for r ak > ^ak(WSC)> r a k(SSC) the model is consistent [i.e., equations (5) or (6) are verified], 
this means that the true critical anisotropy radius (r a k c ) for that model must satisfy the relation 
r ak(NC)< r a k c < ^ a k(WSC)< r a k(SSC). In this case, the limitations on r a k obtained from WSC or 
SSC, are upper bounds on the lower limit r a k c for consistency. Obviously, the situation can be 
more complicated [see equations (38)-(39), and the following discussion]. 

In the next section, after presenting the (71,72) models, the analytical constraints on r a for 
the consistency of the one-component OM anisotropic 7 models are derived, together with a 
general result on the consistency of two-component isotropic (71,72) models. Successively, more 
specific results for the isotropic (1,0) models are proved. Finally, a limitation on r a for the OM 
anisotropic (7,BH) models, is explicitly derived as a function of 7. 



3. The (71,72) Models 



Both density distributions of the (71,72) models belong to the widely explored family of the 
so-called 7 models ( pehnen 1995 , hereafter D93; jCarollo 1993| , Tremaine et al. 1994): 



p(r) 



7 



Mr r 



4vr rT(r c +r) 4 -T' 



M(r) = M x 



r c + r 



3-7 



< 7 < 3, 



(7) 



where M is the total mass and r c a characteristic scale-length. The corresponding relative 
potential is given by 



ij/( r ) 



GM 



»c(2-7) 



1 



r 



r + r c 



2-7' 



*(r) 



GM, r + r c 
In 



(8) 



where the second expression holds for 7 = 2. In the following, the mass M = M 71 and the 
characteristic scale-length r c = r c i of the 71 model will be adopted as normalization constants, so 
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that from equation (7) it follows p 7l (r) = p^p 11 (s) and p l2 {r) = p^pp^ 2 (s, (3), where s = r/r c i, 
Pn = -^7i/ r ci5 A* = M 72 /M 7l , an d r c2 = /S^ci- The fundamental ingredient in recovering the 
DF is the total potential V&t = ^71 + ^72) where from equation (8) ^ 71 (r) = ^n^ti ( s ) an d 
^ 72 (r) = \I , n/U , I' 72 (s, (3), and V&n = GM 71 /r c i. With this choice, the (71,72) models are structurally- 
determined by fixing the four independent parameters (M 7l ,r c i,//,/3), with the obvious condition 
// > and /3 > 0. 

For future reference, I give here the explicit expressions of the density and the 
potential for the (1,0) models, for which (M 71 , r c i , p 71 , * 71 ) = (Mi, ri, pi, \I/i), and 

(M 72 ,r c2 ,p 72 ,^ 72 ) = (M ,r , p ,^o): 

Pi(r)=p N/ 5i( s) = g-^l_, (9) 

and 

Pair) = PnWo(s) = H 4?r TgTj ^ 4 , (10) 
where s = r/r±, p^ = M\/r\, Mq = pMi, and ro = moreover, 

*i(r) = Mi(s) = (11) 
* (r) = ^NP^o(s) = ^ 2( ^ 2 g g )2 , (12) 

with ^ N = GM 1 /r 1 . 



3.1. The Necessary and Sufficient Conditions for the 7 Models 

Here I study first the NC for the general case of the anisotropic one-component 7 models, in 
order to determine analytically a critical anisotropy radius such that a higher degree of radial OM 
anisotropy (i.e., a smaller anisotropy radius) produces a negative DF for some permitted value of 
Q, no matter what kind of halo density distribution is added. The unit mass and unit length are 
the total mass M and the scale-length r c of the model, with s a = r a /r c . 

As shown in Appendix A [equations (A1)-(A2)], for 2 < 7 < 3 the NC is satisfied for s a > 0, 
i.e., the possibility that 7 models with 7 > 2 are assembled using only radial orbits is left open by 
the NC. On the contrary, for < 7 < 2 the NC requires 

/ 2-7-2s M 

* a - gM V 7 + 4*m ' (13) 

where sm = sm(j) is given by equation (A2). In this case the NC proves that 7 models with 
< 7 < 2 cannot sustain radial orbits only. In Fig. 1 the lower bound for the anisotropy radius 
as a function of 7 derived from the NC is shown. From the discussion in Section 2, it follows that 
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all 7 models (one or multi-component) in the nearly triangular region under the solid curve are 
inconsistent. 

The WSC can be treated analytically for the one-component 7 models, as shown in Appendix 
A [equations (A3)-(A6)], and we obtain the following limitats on s a : 



3/2 / j ~ 7 ~ S M ( u) 

~ Sm « 6^ + 2(1 + 7)^ + 7' 



where sm = sm(7) is given by equations (A4)-(A6). As discussed in Section 2, the r.h.s. of 
equation (14) (represented in Fig. 1 by the dotted line), is an upper limit on the lower bound for 
the critical anisotropy radius as a function of 7: all one-component 7 models in the region above 
the dotted line are consistent. 

A stronger limitation on r a is obtained using the SSC, but unfortunately this condition for a 
generic 7 results in a trascendental equation that cannot be solved explicitly. However, as shown 
in Appendix A, for the three values 7 = (0, 1, 3) the solution can be derived explicitly. For 7 = 0, 



3(l + 2s M -s 2 
Us^ + Ws M + 2 



fla> W T7V Tn " M o -0-501, (15) 



where sm = sm(0) is given by equation (A8). For 7 = 1, equation (A9) shows that 



s a >4\ -^_ 0.250, (16) 
y 28s^ + 17sm + 4 

where sm = sm(1) is given by equations (AlO)-(All). The numerical application of the SSC to 
the 7 = 2 model flJaffe 1983| ) gives Sa^0.047. Finally, the case 7 — 3 is trivial, the SSC reduces to 



s& > 0. These four values are represented in Fig. 1 by black dots: all one-component 7 models in 
the region above the dashed line are consistent. Note how the SSC improves the estimate of the 
lower bound of r a with respect to the WSC. 

As already pointed out in Section 2, the true critical value of s a for any one-component 7 
model is between the NC and the SSC curves. For the 7 = and 7=1 models these values, 
determined directly from their DFs (Sections B1-B2 in Appendix B), are ~ 0.445 and ~ 0.202, 
respectively. Merritt (1985b) derived the analytical DF for a totally radial Jaffe model (the 7 = 2 
case): its positivity implies that in this case the true lower limit on s a is zero. In Fig. 1 these three 
values are represented by black squares, and in Table 1 all the previous results for the specific cases 
7 = (0, 1,2,3) are summarized: note how the more the model is concentrated, the more radial 
anisotropy can be supported. The true critical r a value for the OM anisotropic one-component 7 



models as a function of 7 is numerically known (see Fig. 1 in Carollo et al. 1995a ). 
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3.2. Sufficient Conditions for Isotropic (71,72) Models 

In order to proceed further with this analytical discussion, and to allow for the presence of a 
"halo" component, we will use the WSC rather than the more complicated SSC. The following 
three results are proven analytically in this Section: 

1. In the case of globally isotropic two-component (71,72) models with 1 < 71 < 3 and 

< 72 < 71, the DF of the more peaked component 71 is positive over all the phase space, 
for all values of (fi, (3) = (M 7l /M 72 , r C 2/r c i)- As a consequence, the 7 = 1 component of (1,0) 
models is consistent for all values of the parameters (/z, (3). 

2. In the case of globally isotropic (1,0) models, the WSC applied to the 7 = density 
distribution suggests the existence of a lower limit of fj, = Mq/M\ as a function of /3 = tq/tx. 
In particular, for (3 < 5/2, i.e., when the 7 = component is sufficiently concentrated, all 
values of ii can be accepted. Using the analytical DF the existence of this lower limit will be 
proved in Section 5. 

3. In the case of anisotropic 7 models with a BH at their center, it is possible to determine 
analytically a lower limit on r a as a function of 7 using the WSC. 

The proof of the first result is conceptually straightforward but algebraically cumbersome. 
In Appendix A [equations (A12)-(A13)] it is proved that under the hypothesis assumed in point 
1 above, equation (6) is verified for all choices of (/i,/5). I note explicitly that this result contains 
as a particular case the fact - already proved in C96 - that globally isotropic (1,1) models can be 
self-consistently assembled for any choice of (/U,/3). From this general result it also follows that - 
e.g. - the same is true for isotropic Jaffe+Jaffe models (71 = 72 = 2). Finally, considering that for 
r c — > the potential of 7 models becomes that of a point mass [see equation (8)], the previous 
result means that a BH of any mass can be added at the center of a globally isotropic 7 model 
when 1 < 7 < 3. A different analysis, based on a series expansion of the integral representation of 
the DF for isotropic 7 models, shows that the true limit on 7 in order to allow for the presence of 
a BH of any mass at the models center, is 7 > 1/2 ( [Tremaine et al. 1994| ). 

The result stated in point 2 above is proved in Appendix A [equation (A14)]. For equation 
(6) to be satisfied for the isotropic 7 = component of (1,0) models, a sufficient condition for the 
consistency of this component is 

fi> (2/3-5)/3 2 . (17) 

This requirement can be interpreted in two different ways. The first is that, having fixed the ratio 
[3 = vq/ti, only for sufficiently high mass ratios /x = Mq/M\ can the 7 = component "dilute" 
the effect of the central cusp of the 7 = 1 model on the total potential, and be consistent. More 
specifically, when /3 < 5/2 (i.e., when the 7 = density distribution is sufficiently concentrated), 
even a vanishing mass Mq (fj, — > 0) can be accepted, while for large (3 only very large \x are allowed. 
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From another point of view, equation (17) tells us that having fixed fx, (3 cannot be arbitrarily 
large, but in some sense the concentration of the 7 = component must adapt to the density 
distribution of the 7=1 component. The effect of the concentration is much more important than 
the amount of mass: in fact /3^(/^/2) 1//3 . This means that even increasing considerably the mass 
ratio, the maximum value of ro allowed for the 7 = model grows only like the cube root of /x. 
The limitation j3 < 5/2 is only a sufficient condition for the consistency of a 7 = model coupled 
with a dominant 7 = 1: a larger critical value for (3 is expected from direct inspection of the DF 
when \i — > (see Section 5). 

The result presented in point 3 above can be interpreted as an extension to the radially 
anisotropic case of the analysis performed by Tremaine et al. (1994), and is proved in Appendix 
A by showing that the WSC applied to the anisotropic (7,BH) models with 1 < 7 < 3 can 
be analytically discussed in the special case of a dominant BH, i.e., assuming in equation (6) 
Mp = Mbh (and so = GMbk/t). Unfortunately, in the non-asymptotic case, the equation 
to be discussed is transcendental, and no analytical discussion can be carried out. At first, the 
assumption of a dominant BH could appear as a very rough approximation of reality, but this 
is not true: the constraint derived can be used as a safe limitation when constructing models 
containing a BH of a realistic mass at their centre. As shown in Appendix A, [equation (A15)], for 
1 < 7 < 3, 

, > , / (3-7)(7-2)+4(3-7)*M-247 

Sa " SM V 12^ + 8( 7 -1)*m+ 7(7-D ' (18) 
where sm = «m(7) is obtained by solving a fourth degree algebraic equation. In Fig. 1, the 
long-dashed line represents the lower bound for s a as determined by the previous equation, while 
the explicit values for 7 = (1,2,3) are given in Table 1. In particular, the critical value for the 
7 = 1 model was already derived as a limiting case of a (non-asymptotic) formula given in C96 
[equation (15) there]. 



4. The DF of (1,0) Models 

We can now proceed to the explicit recovery of the DF of the (1,0) models. Just as for 
the density and the potential, it is also useful for the DF to work with dimensionless functions; 
the two components of the DF are of the form / = /n/O^AQ) with /n = p^^^ 3 ^ 2 and 
< Q = Q/^n < ^i(0) + /U^o(0). The easiest way to compute each DF is to use the first of 
the identities in equation (2). For the evaluation of the integral one would be tempted to obtain 
g(^f^), eliminating the radius from the modified density and the total potential: formally, this can 
be done, but the resulting expression for the radial coordinate involves a quadratic irrationality, 
that, after insertion in equations (3), (9) and (10), produces an intractable expression. Here I 
follow another approach: instead of eliminating the radius, the integration variable is changed 
from the total potential to the radius itself. This is equivalent to a remapping of the domain of 
definition of each /, from the range of variation of to the range of variation of r, and leads us 
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Q = -T—.+v Jta : ,A2 ' 0<^<CX). (19) 



to introduce the dimensionless radius v using equations (11)-(12): 

1 (3 + 2v 

As shown in Appendix B [equations (B1)-(B2)], with this change of variable and after normalization 
to the dimensional scales of the Hernquist density distribution, the DF for the 7 = 1 and 7 = 
components can be formally written as: 



f(Q) = MQ) + UQ) - /x 




si 



u = u(Q), (20) 



where s a = r a /r±, and the subscripts refer to the isotropic and anisotropic parts of the DF, 
respectively. Following the procedure, f(Q) results from the elimination of v between equations 
(19)-(20). In the general case, i.e., for any choice of (fi,/3, s a ), f(Q) [and the so called differential 
energy distribution for each component as well, (see, e.g., BT87, p. 242)] can be recovered 
analytically [see equation (B3) for a proof of this fact]. So, it is shown that in addition to the (1,1) 
models, the (1,0) models are a class of two-component stellar systems in which both the spatial 
density distributions, the solution of the Jeans equations (see Appendix C), and the phase-space 
distribution functions can be explicitly found. Unfortunately, the DF of the (1,0) models results in 
a combination of elliptic functions, even more complicated than the DF of the (1,1) models, and 
this limits their applicability to special problems in which the DF is required to be known with 
arbitrary precision or to be formally manipulated. 

Here I present only the DFs for the two density distributions obtained under the assumption 
of a dominant "halo" component; I derive the DF for a 7 = 1 model with a dominant 7 = halo 
(// — > 00), and for a 7 = model with a dominant 7 = 1 halo (/x — > 0). Technically, this reduces 
to the assumption that the total potential is the potential of the halo component only. Even 
though it is a limiting case, the study of halo -dominated models is interesting for several different 
reasons: 1) the formulae - expressible using elementary functions - are much simpler than in the 
general case, and can be studied very easily, making clearer the effect of the halo component on 
the DF; 2) the halo-dominated case is the one that differs most from the case of the corresponding 
one-component model, and so the differences are better evident; 3) all the intermediate cases 
fall between the one-component model and the halo-dominated one. A comparison with more 
realistic values of the halo masses is postponed to Section 5. 

In the following paragraphs, the two DFs will be compared with those of the corresponding 
one-component 7 = 1 and 7 = models, and the exact phase-space constraints will be derived 
and compared with those obtained using the NC, WSC, and SSC in Section 3. 



4.1. The 7 = 1 Model Plus a 7 = Dominant Halo 



The explicit expression for the DF of the 7 = 1 model with an arbitrary degree of OM orbital 
anisotropy immersed in a dominant 7 = halo is derived here. Formally, this case corresponds to 
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the assumption of fi — > oo in the total potential, i.e., \£t = ^o, and so in equations (19)-(20) 

After differentiation inside the integral in equation (B2) with g given by equations (3) and (9), 
and after a partial fraction decomposition of the rational part of the integrand, one obtains: 

F[(v) = ^±^==\Rl + ml - Hi - (1 + ml - (1 + W)H\ - 3(0 - l)H\], (22) 

and 

* M = ^m {ml - (5 - 2,ml - 3W - 1)Ht] - (23> 

The H functions depend on [3 and u, and are defined as 



£ = ^, and X = -Ar. (25) 



where 



When £ = 1 and n > 1 

zrz m _ v^r(n) 

" r(n-l/2)(«/ + A)»' (26) 
where T is the complete gamma function. When £ 7^ 1 the recursion formula 

ffZ (n = 1 dHZiQ _ (-1)" cPfff(g) 

n(^ + A) d£ n!(i/ + A) n d£» 1 j 

holds, and so the explicit evaluation of Hf(£) suffices: 

2 [^^' ifO<£<l; 

iVS/ IV + A I arccoshy^ f £ V ' 

In order to distinguish between the two cases £ > 1 and < £ < 1, a careful discussion is needed. 
From equation (25) the value z = corresponds to £ = ^/(f + A) < 1 V(i>,/3), and so the first 
of equations (28) must be used for the evaluation of all H® functions. More complicated is the 
case z = 1, when £ = (1/ + + A): note that for (3 > 0, A is a monotonically increasing 

function of v, with A = for v = and A — > /3/2 for 1/ — > 00. As a consequence, it follows that 
W, < /3 < 2 £ > 1. When instead, > 2, 3 z/ cr = /?/(/? - 2) so that 1/ < i/ CT £ > 1, 
f = v cv £ = 1, and v > u cv £ < 1. This completes the derivation of the DF for the 7 = 1 
component. 

In Fig. 2 the comparison with the DF of the one-component 7=1 model (solid line), in case 
of global isotropy and for a specific value of the anisotropy radius, is given. Such DF was given in 
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H90 as a function of Q, but for consistency with the present work it is derived in Appendix B as a 
function of v [equations (B4)-(B6)]. The formulae derived in this paragraph have been tested for 
many values of (3 and r a using a code that numerically recovers the DF for spherically symmetric 
multi-component galaxy models with OM anisotropy; we obtained extremely good agreement; in 
all cases the maximum differences between the analytical and numerical DFs are much less than 1 
per cent. 

In the upper panel of Fig. 2 the isotropic case is presented. Note how for (3 > 1 the DF 
is more peaked than for the one-component 7 = 1 model, and the opposite holds when (3 < 1: 
this behaviour was already found in the isotropic (1,1) models (C96, Fig. 2). In the lower panel 
the anisotropic case is shown when s a = 0.26, near the consistency limit for the one-component 
7 = 1 model (see §3.1). The main effect of anisotropy, as already found for (1,1) models and 
Ri/rn models (CL97), is the appearance in the DF of a depression well outside the galaxy center. 
Decreasing the anisotropy radius, the depression deepens, running finally into negative values for a 
critical value of s a (dependent on [3) and making the model inconsistent. Again, as already found 
for (1,1) models, this effect is stronger for smaller (3 values, i.e., a very concentrated halo makes 
the DF more sensitive to the effects of anisotropy, while the opposite is true for halos more diffuse 
than the 7 = 1 density distribution. 



4.2. The 7 = Model Plus a 7 = 1 Dominant Halo 

The explicit expression for the DF of a 7 = model with an arbitrary degree of OM orbital 
anisotropy, immersed in a dominant Hernquist halo, is derived here. Formally this case corresponds 
to the assumption of fi — > in the total potential, i.e., = ^1, and so in equations (19)-(20) 

* = IT7 (l)"-' 1 + "' 2 < 29 > 

After differentiation inside the integral in equation (B2) with g given by equations (3) and (10), 
and after a partial fraction decomposition of the rational part of the integrand, one obtains: 



and 



where 



F i{ u) = 3 J^±^4G 5 , (30) 

47T 

Am = 3l " 3 ^ T ^ (tf 2 °5 - wg, + 2G 3 ), (si) 



When = 1, and n > 2 

V^Tjn- 1) 
r(n-l/2)(l + i/) 



Gn{i^) = — v-y;- ti__ . (33) 
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When (3^1 and n > 2 the recursion formula 



G n+ i{p,v) 



I dGJJ^v) = j-iy- 1 d n ~ l G 2 {P,v) 
n dp n\ dj3 n ~ 1 



(34) 



holds, and so the explicit evaluation of G 2 suffices: 



P + v 



1 



+ 



(P + vfl 2 



1 + v 




if < < 1; 



if > 1. 



(35) 



In Fig. 3 the comparison with the DF of the one-component 7 = model (solid line), in the case 
of global isotropy and for a specific value of the anisotropy radius, is given. Such DF was given 
in D93 as a function of Q, but for consistency with the present work it is derived in Appendix B 
as a function of v [equations (B7)-(B9)]. As in the previous case, the derived formulae have been 
successfully tested for many values of P and r a by comparison with the numerically derived DFs, 
obtaining maximum differences less than 1 per cent in any case. 

In the upper panel the isotropic case is presented. Note how for P < 1 the DF is more peaked 
than for the one-component 7 = model, and the opposite holds when P < 10. This behavior 
is similar to that found in the previous section. In the lower panel the anisotropic case is shown 
when r a = 0.65ro, near the consistency limit for the one-component 7 = model (see §3.1). As for 
the 7 = 1 component, the main effect of anisotropy is the appearance in the DF of a depression 
well outside the galaxy center, and again the depression becomes deeper and deeper decreasing 
the anisotropy radius. Finally, as in the previous case, the effect of anisotropy is found to be more 
important when the halo concentration increases. 



We move now to comment on the main similarities and differences between the DFs of the 
two components of the (1,0) models, especially considering the role of concentration and orbital 
anisotropy in determining their consistency. For simplicity the discussion is restricted to the 
halo-dominated cases. 

The first important point addressed by using the DFs, is the study of the effect of halo 
concentration in determining the consistency of the two (1,0) model components in the case of 
global isotropy. The effect of the 7 = halo concentration on the consistency of the globally 
isotropic 7 = 1 component, can be derived by the direct inspection of the DF and confirms 
the analytical prediction obtained using the WSC in §3.1. In fact, it is found that the globally 
isotropic 7 = 1 component is consistent independent of the concentration and total mass of the 



2 Having defined (3 = ro/n, at variance with what happened in §4.1, a more diffuse 7=1 component corresponds 
to /3 < 1, and vice-versa. 



5. Consistency of (1,0) Models 
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superimposed 7 = halo: only anisotropic 7 = 1 component in (1,0) models can be unphysical 
due to the presence of the 7 = density distribution. 

For the globally isotropic 7 = component with a dominant 7 = 1 halo, the situation is 
more complicate, because, in accordance with the analysis presented in §3.2, its DF may become 
negative in case of a high concentration of the external 7 = 1 component. In fact, the DF becomes 
negative for £ — > x I / i(0) when /3>5.233, a larger value than the more conservative one (5/2) 
derived using the WSC. A closer look at this behavior, and a comparison with the qualitatively 
different behavior exhibited by the DF of the 7 = 1 component, is particularly instructive. In 
fact, while the DF of the 7 = 1 density distribution diverges at high (relative) energies both 
in the one-component and in the halo-dominated cases (Fig. 2), the DF of the 7 = model 
is divergent for high energies in the one-component case, but finite in the halo-dominated one 
(Fig. 3). Moreover, when increasing the 7 = 1 halo concentration (i.e., increasing = ro/77), the 
central value of the DF associated with the 7 = density profile decreases monotonically, and, 
for greater than the before mentioned critical value, it becomes negative, revealing the model 
inconsistency. It must be stressed that a similar behavior was found in the numerical investigation 
of consistency of King (1972) and quasi-isothermal halos added to a de Vaucouleurs (1948) density 
distribution, carried out by CP92. Also, note how the decrease of the central value of the DF for 
increasing halo concentration is reminiscent of that found by C96 for (1,1) models, even if in that 
case the transition was found to be more discontinuous: the DF of a 7 = 1 component in (1,1) 
models remains divergent at the center for all finite concentrations of the other 7 = 1 component, 
and becomes exactly zero at the center only when the halo is reduced to a central BH (see Fig. 2 
in C96). The qualitative discussion above can be put on more quantitative grounds. In fact, in 
the halo-dominated case, the central value of the DF of the 7 = component is easily derived for 
a generic using the formulae given in §4.2: 



3/U 



# = _^_<f, if = 1, (36) 



3(64/3 3 -240/3 2 +280/3-105) „ rpt „ n /T^3 16 / 3 3 -328 / 3 2 +630/3-315 - f n R -, 
32(1-/3)5/2/39/2 arctany — 3 2(i-/3) 2 /3 4 > u u < p < i, 



5 ' 

3(64/? 3 -240/3 2 +280/3-105) , i /Tj^T 16/3 3 -328ff 2 +6300-315 - f R , . 
32(/3-l) 5 /2/3 !) /2 dICtdliny ^ 32(/3-l) 2 /3 4 ' 11 P > X > 



and 



J a 



3/x 



8^/271-3 



3(8/3 2 -12/3+5) 



if < < 1, 



32(1-/3) 

f, if = 1, (37) 

3(8/3 2 -12/3+5) arrfan u //EI + (4/3-3) (2/3-5) . f g -, 
, 32(/3-i) 5 /2 /3 5/2 arctanny ^ -f ^ 2 {p-\) 2 2 > 11 P > i - 

The limiting value for [3 in the isotropic case is obtained by solving numerically the equation 
/° = for > 1. In Fig. 4 (where the high-concentration case corresponding to > 1 is shown), 
the decrease of /° when increases is apparent. 

The second important point addressed by using the DFs, is the study of the combined effect 
of orbital anisotropy and halo concentration in determining the consistency of the two (1,0) model 
components. We cannot expect a simple behavior, because - as should be clear from the previous 
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sections - halo concentration and anisotropy affect the DF in different regions of the phase-space, 
i.e., the high energy regions of the DF are more sensitive to concentration effects, while the OM 
orbital anisotropy acts mainly at intermediate energies. The simplest way to summarize the results 
is to express the consistency limitations in terms of the anisotropy radius of each component as 
function of /3, determining in the parameter space (s a , (3) the critical regions where the models 
are consistent. This approach is particularly useful because independent of the specific form of 
the density profile of the investigated model, the positivity requirement for each DF of an OM 
multi-component system over all the phase-space can be expressed in term of the anisotropy 
radius as a function of the other model parameters, due to the simple appearance of s a in equation 
(20). 

In fact, let be A + the set defined by the property that /; > W G A + . Then, from equation 

(20), 





max |o, sup 


[ /aM" 

















(38) 

is a first condition to be satisfied. Obviously, when /; > over all the phase-space (the common 
situation), A + coincides with the total range of variation for u, and equation (38) is also the only 
condition to be checked for the model consistency. In this case equation (38) shows that there is 
at most a lower bound for the anisotropy radius, s~ . For example, this is the case for the 7 = 1 
component in (1,0) and (1,1) models, or for one-component anisotropic 7 models. 

When the set (complementary to A + ) is not empty, i.e., /; < over some region of 
phase-space, a second inequality, derived from equation (20), must necessarily be verified: 



5a < Sac = 



\ 



inf 



(39) 



A general consequence of equations (38)- (39) valid for all single or multi-component spherically 
symmetric, radially anisotropic OM models, is that the allowed region for consistency in the 
anisotropy space is given by s~ < s a < s+. Moreover, if / a < W € A-, or s+ < s ac , then the 
proposed model is inconsistent. 

The quantitative trend of s~ c for the 7=1 density distribution with a dominant 7 = halo 
is shown in Fig. 5a (solid line): for a given (3 all values of s a higher than the critical curve are 
acceptable. Note how an increase in the halo concentration (a decreasing (3) produces an increase 
of s ac , i.e., a very concentrated halo makes the other component more sensitive to anisotropy 
effects, a behavior qualitatively anticipated in §4.1, and already found for (1,1) models (C96, Fig. 
5). 

A more complicated (and more interesting) case is presented by the halo-dominated 7 = 
model. In this case we already know that, due to the halo concentration, even the isotropic 
case can be inconsistent, i.e., /; < 0. This means that s ac must also be considered. The trend 
of s~ c is shown in Fig. 5ab (dotted line): as in the previous case an increase of the minimum 
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anisotropy radius corresponds to an increasing halo concentration (i.e., to an increasing j3). As 
(3 increases above the critical value ~ 5.233, fi becomes negative at the center, and the isotropic 
7 = component becomes inconsistent: in Fig. 5a this region is contained in the box at the 
top-right, that is enlarged in Fig. 5b. Here the dotted line is again s~ c , and the dashed line 
represents s+: for 5.233^/3^6.15 the inequality s~ c < s+ holds, and according to equations 
(38)-(39) the region between the two curves corresponds to consistent 7 = components. This is 
a quite counterintuitive example of the combined effect of an external potential and anisotropy on 
the consistency of an anisotropic galaxy model, where an otherwise inconsistent isotropic model 
is made consistent by orbital anisotropy! Finally, for /3>6.15 no physically acceptable 7 = 
components are possible, even considering the effect of anisotropy. 

A question arises: how well does the asymptotic analysis obtained in the limit of dominant 
halos compare to the more realistic cases of halos with finite mass? An answer can be obtained 
by inspection of Fig. 5a, where the dashed lines represent the limits on the anisotropy radius 
obtained when considering a halo ten times as massive as the component investigated. Note that 
when the halo is more concentrated than the considered density component [large (3 in case of 
7 = model (dotted line) and small (3 for the 7 = 1 model (solid line)], the curves corresponding 
to the asymptotic analysis and the dashed ones are indistinguishable for any practical application. 
On the contrary, a small departure appears when the halo scale-length is substantially larger than 
that of the considered density component, with the dashed curves approaching the critical value 
for the anisotropy radius corresponding to the one-component model (the two black dots). This 
is an obvious behavior, since for any finite value of the halo mass, its gravitational effect becomes 
weaker and weaker for larger and larger halo scale-length. 

6. Conclusions 

In this paper, an extensive analytical investigation of the phase-space of two-component 
spherical galaxy models made of the sum of a Hernquist density distribution and a 7 = 
model with different physical scales, is carried out. Following the simple Osipkov-Merritt 
parameterization, a variable amount of orbital anisotropy is allowed in each component. For these 
models, other important properties useful in applications - the velocity dispersion components 
and the various energy terms entering the scalar virial theorem - can be expressed analytically, 
and are given in Appendix C. The main results can be summarized as follows: 

1. The necessary and sufficient conditions that the model parameters must satisfy, in order 
to correspond to a (1,0) system for which the two physically distinct components have a 
positive DF are analytically derived using the method introduced in CP92. Some conditions 
are obtained for the wider class of two-component (71,72) models [of which the (1,0) models 
are a special case]. In particular, it is shown that the DF of the 71 component in isotropic 
(7i)72) models is nowhere negative, independent of the mass and concentration of the 72 



-18- 



component, whenever 1 < 71 < 3 and < 72 < 71. As a special application of this result, 
it follows that a BH of any mass can be consistently added at the center of any isotropic 
member of the 7 family of models, when 1 < 7 < 3. Two important consequences follow. The 
first is that the consistency of isotropic (1,1) [or (1,BH)] models proved in C96 using an "ad 
hoc" technique is not exceptional, but a common property of a large class of two-component 
7 models: for example, also isotropic two-component Jaffe (7 = 2) or Jaffe+BH models can 
be safely assembled. The second is that in two-component isotropic models, the component 
with the steeper central density distribution is usually the most robust against inconsistency. 

2. It is shown that an analytic estimate of a minimum value of s a for one-component 7 models 
with a massive (dominant) BH at their center can be explicitly found. As expected, this 
minimum value decreases for increasing 7. 

3. It is shown that the analytic expression for the DF of (1,0) models with general OM 
anisotropy can be found in terms of elliptic functions. The special cases in which each one of 
the two density components are embedded in a dominant halo are also discussed: under this 
assumption the DFs can be expressed using just elementary functions, allowing a detailed 
analytical investigation. 

4. The region of the parameter space in which (1,0) models are consistent is explored using the 
derived DFs: it is shown that, unlike the 7 = 1 component, the 7 = component becomes 
inconsistent when the halo is sufficiently concentrated, even in the isotropic case. This is 
an explicit example (albeit not so extreme) of the result found by CP92, that numerically 
proved the impossibility of adding a King or a quasi-isothermal halo to a de Vaucouleurs 
galaxy. In such models, the (isotropic) de Vaucouleurs galaxy was found instead consistent 
over all the parameter space. 

5. The combined effect of halo concentration and orbital anisotropy is finally investigated. The 
trend of the minimum value for the anisotropy radius as a function of the halo concentration 
is qualitatively similar in both components, and to that found for (1,1) models in C96: a 
more diffuse halo allows a larger amount of anisotropy. A qualitatively new behavior is found 
and explained by investigating the DF of the 7 = component in the halo-dominated case 
for high halo concentrations. It is analytically shown that there exists a small region in the 
parameter space where a sufficient amount of anisotropy can compensate the inconsistency 
produced by the halo concentration on the structurally analogous - but isotropic - case. 

As a final remark, it can be useful to point out some general trends that emerge when 
comparing different one and two-component models with OM anisotropy, as those investigated 
numerically in CP92 and CL97, and analytically in C96 and in this paper. The first common 
trend is that OM anisotropy produces a negative DF outside the galaxy center, while the halo 
concentration affects mainly the DF at high (relative) energies. The second is that the possibility 
of sustaining a strong degree of anisotropy is weakened by the presence of a very concentrated 
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halo. The third is that in two-component models, in cases of very different density profiles in the 
central regions, the component with the natter density is the most "delicate" and can easily found 
to be inconsistent: particular attention should be paid in constructing such models. 

I would like to thank Giuseppe Bertin, Laura Greggio, and Silvia Pellegrini for helpful 
comments and discussions. The referee, Stephen Levine, is especially thanked for his comments that 
greatly improved the paper. This work has been partially supported by contracts ASI-95-RS-152, 
ASI-ARS-96-70, and MURST-Cofin98. 



A. Consistency Requirements 

A.l. The NC and WSC for One Component 7 Models 

The NC [equation (4)] for the anisotropic (one-component) 7 models imposes a limitation on 
the anisotropy radius: 

2 s 2 (2-7-2s) IK . 

si > — -, -, < s < 00. (Al) 

a ~ 7 + 4s ' ~ ~ K J 

This inequality is true for s 2 larger than or equal to the maximum of the function on the r.h.s. 

of equation (Al). For 2 < 7 the r.h.s. is strictly negative, and so all values of s a satisfy the NC. 

When < 7 < 2 the maximum is reached at 



s m(7) ' ' > 

that after substitution in equation (Al), gives equation (13). 

The WSC [equation (6)] applied to the anisotropic one-component 7 models gives the 
following inequality: 

9 s 3 (3 — 7 — s) ... 

s 2 > „ v . — '—— - , < s < 00. (A3) 

a ~ 6s 2 + 2(l + 7 )s + 7' - - v ; 

After differentiation, one is left with the discussion of a cubic equation. Its discriminant is negative 

for < 7 < (-\/73 - 5)/8 and positive for (\/73 - 5)/8 < 7 < 3. In the first case two of the three 

real solutions are negative, and the maximum of the r.h.s. of equation (A3) is reached at 



1-7 2^4^ 
sm{7) = — — + cos 



1 yp=w=w=SF*?J 



3 II + II7-47 2 

In the second case, discarding the two complex conjugates roots, and defining 



(A4) 



s = (4 - 7 )[11 + 11 7 - 4 7 2 + yj(15 - 47) (3 - 7 )(47 2 + 57 - 3)], (A5) 
the maximum is reached at 

SM(7) .£f! + I_l + ^. (A6) 



Finally, when 7 = (V73- 5)/8 ~ 0.443, s M = (V73 + 3)/8 ~ 1.443 
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A.2. The SSC for the 7 = and 7 = 1 Models 

The SSC [equation (5)] applied to the anisotropic one-component 7 = model gives the 
following inequality: 

2 3s 2 (l + 2s-s 2 ) , A . 

sl > — h -, < s < 00. (A7) 

a ~ 14s 2 + 10s + 2 ' ~ ~ v ; 

The maximum of the r.h.s. of equation (A7) can be obtained explicitly solving a fourth-degree 
equation. The numerical value of the only physically acceptable solution is 

s M (0) ~ 1.3149, (A8) 

that after back substitution in equation (A7) gives equation (15). 

The SSC applied to the anisotropic one-component 7=1 model gives the following inequality: 

2 3s 3 (3-2s) /A x 

sl > ^ — z — ~~r, < s < 00. (A9) 

a - 28s 2 + 17s + 4' ~ ~ v ; 

After differentiation, discarding the two complex conjugates roots of the resulting cubic equation, 
and defining 

s = 681939 + 84^35887965, (A10) 
the maximum of the r.h.s. of equation (A9) is reached at 

s 1 / 3 3 1987 

«« = w-| + ^ = a9116 - < A11 > 

that after back substitution in equation (A9) gives equation (16). 



A. 3. The WSC for Isotropic (71,72) 

With the aid of the WSC, the positivity of the DF for the 71 density distribution of the 
globally isotropic two-component (71,72) models, where 1 < 71 < 3 and < 72 < 71, is here 
discussed. Having computed the derivatives prescribed by equation (6), we have to investigate the 
positivity of a rational expression, whose denominator is strictly positive V(7i,72) and V(s,/3); the 
numerator factorizes in a strictly positive function and in a transcendental expression. Defining 
71 = 1 + ei (with < e\ < 2) and 72 = 71 — £2 (with < 62 < 71), the trascendental factor reads: 

2(s + f3)\l + 1/s) 71 [6s 2 + 2(1 + 7 i)s + 71] + /i(a + 1) 3 (1 + f3/s)^F(s, (3, 71, 72), (A12) 

where the first addend is strictly positive, and 



F(s, 13, 71 , 72) = 12s 3 + 4[(5 - ei + e 2 )P + 2ei]s 2 + [(10 + 5ei - e\ + 5e 2 + e^P + em]* + /3 7 i (2 + e 2 ) . 

(A13) 
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In the range of values for e\ and 62 the positivity of F is easy proved V(s, /3) > 0, and so condition 
(6) is verified V(/t,/?). 

The application of the WSC to the globally isotropic 7 = component of a (1,0) model is 
simple, leading to the discussion of: 

-(3s + 5-2/3)(s + /3) 3 , K 

" s (3. • - s -°°- (A14) 

First, note that for [5 < 5/2 the previous inequality is satisfied V/i > 0. When (5 > 5/2, the 
maximum of the r.h.s. of equation (A14) over the interval < s < 00 is reached at s = 0, and so 
the condition given in equation (17) is derived. 



A.4. The WSC for Anisotropic ( 7 ,BH) Models 

In this case we assume in equation (6) Mp = Mbh, and, from the previous discussion, 
we restrict to 7 models with 1 < 7 < 3. After the computation of the derivatives, we have to 
investigate the following inequality 

^[(3--,)( 7 -2) + 4(3- T )»-2^] 

12s 2 + 8(7- 1)5 + 7(7- 1) ' - - v 1 

After differentiation, one is left with the discussion of a quartic equation, and it can be shown that 
there exists only one maximum, located at sm = sm(i) > 0. The explicit expression for sm(t) is 
not very useful, and so is not reported here. In any case, it can be of interest the explicit values of 
sm(7) for 7 = (1,2,3). After easy computations, one finds 

s M (l) = 2, s M (2) = (M + 6 f* )1/3 + — - 6 2 ^_ )i/3 ~ 1.191, SM (3) = 0, (A16) 

that after substitution in equation (18) gives the values reported in Table 1. 



B. The DF of (1,0) Models 

Here the main steps required for the determination of the DF of each component of (1,0) 
models are briefly described. Changing the integration variable from the total potential to the 
radius, and after normalization to the physical scales of the 7 = 1 component, equation (2) 
becomes: 

where the relation between v and Q is given by equation (19), and 

pw=-rf -r- ds - (B2) 
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Note that a negative sign appears in front of the previous integral, due to the monotonic decrease 
of the relative potential with radius. From equations (11)-(12) 

i /iT7 v / 2(r+^(/3 + ^)(/3 + g ) 

y/* T (v) - * T (s) V S ~ V VAs^ + Bs + C 

where A,B,C are constants depending on (//,/?, v). Finally, observing that for each component 
dg/ds is a rational function of s, with the standard change of the integration variable 
(1 + s)/(s — v) = t 2 , it follows that the integral in equation (B2) can be expressed in terms of 
elliptic integrals. 



B.l. The DF of the One Component 7 = 1 Model as a Function of v 

In this case, in equation (Bl) 

1 



Q 



1 + 1/' 



dQ 
dv 



(B4) 



After differentiation under the integral in equation (B2), with g given by equations (3) and (7), 
and = ^1, and after a partial fraction decomposition of the rational part of the integrand, one 
obtains: 

2 15^ 2 + 50^ + 59 1 l + i> 



2tt 



15 



(1 + vf 



1 

. -arctan— = 



and 



y/TTv 8 l-5u 
2vr 15 (1 + v) 3 



(B5) 



(B6) 



A check of the derived formulae is obtained by substitution of v = v{Q) as given by inverting the 
first of equations (B4), and comparing the result with the DF given in H90. 



B.2. The DF for the One Component 7 = Model as a Function of v 

In this case, the scale-length ro of the 7 = component and its mass Mq, are assumed as 
physical scales, and so /n = (Mo/ro)(GMo/Vo)~ 3 / 2 . This is formally equivalent to assume (3 = 1 
and jj, = 1 in equations (10) and (12), and so in equation (Bl): 

n_ 1 + 2^ , (dQY 1 _ (1 + ^) 3 (m) 
g -27T+Tp [to) ~ ~- (B7) 

After differentiation under the integral in equation (B2) with g given by equations (3) and (9), 
and = ^0) and after a partial fraction decomposition of the rational part of the integrand, one 
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obtains: 

Fit?) 

and 



3^2 
4tt 



2V1 + 2z/(15z/ 2 + 22z/ + 11) 4(5i/ 2 + 4z/ + 2) 



3\/2 

4"7T 



3(1 + ^) 3 (1 + ^ 



VI + 2v(9v 2 + 2u + 1) 6^ 2 
3(1 + 1/) 3 (l + i/) 



arctanh 



1 



r arctanh 



VI + 2^ 
1 



VTT2^7 



(B8) 



(B9) 



A check of the derived formulae is obtained by substitution of v = v(Q) as given by inverting the 
first of equations (B7), and comparing the result with the DF given in D93fJ 



C. The Velocity Dispersions and Virial Quantities 

Here I present the main dynamical quantities of the model discussed in the paper. The radial 
component a 2 of velocity dispersion a 2 = a 2 + a 2 in the OM parameterization can be written for 
each component as 

ArVAr) =m+i^i (C1) 

where 

A(r) = G r p{r)M T {r)dr, I(r) = G f°° P ^ M ^ dr, (C2) 

J r J r 1" 

( Binney and Mamon 1982| ); once obtained of, the tangential velocity dispersion is given by 

2r 2 

a 2 (r) = ^rf^fAr). (C3) 

Other quantities of interest in applications as the global energies entering the scalar virial theorem 
are derived in the next paragraphs. 



C.l. The 7 = 1 Component 

Due to the presence of the 7 = component, I\ = In + /iiio and A\ = An + fJ,Aio, where 
the dimensional coefficients of the velocity dispersions and energies are \I/n and Un = Mi^Nj 
respectively. After normalization and integration one finds 

~ 1 l + s 12s 3 + 42s 2 + 52s + 25 ~ 1 + 4s 

hl -2^ ln ~s 24tt(1 + S )4 ' Al1 ~ 24vr(l + s)4- (C4) 

The interaction with the 7 = halo is described by the two contributions 

= 3 (3 + s (2s + + l)[6s 2 + 6s(/3 + 1) - (i 2 + 8/3 - 1] 

10 ~ tt(/3 - l) 5 1 + s 4tt(/3- 1) 4 (1 + s) 2 (/3 + s) 2 ' { ' 



3 The anisotropic part of the DF given in D93 should be multiplied by 2. 
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(C6) 



and 

~ (/3 2 + 4/3 + l) 1 (3 + s 

Al0 ~ 2^-1)5 TT7~ 

2{(3 2 + 4(3 + l)s 3 + 3(/3 + l)((3 2 + 4(3+ l)s 2 + 2/?(5/? 2 + 8(3 + 5)s + 6/? 2 (/3 + 1) 

4tt(/3- l) 4 (/3 + s) 2 (l + s) 2 
When /3 = 1 

f 1 I _10s 2 + 5s + l 

/l0 " lMl + a) 5 ' 10 " 60vr(l + • (C7) 

The gravitational energy of the 7 = 1 component is given by U\ = Un + /j,Uio, where 
Un = —2-ir J p\^\r 2 dr is the contribution due to the self-interaction, and U\q = —4ir J p\^or 2 dr, 
is due to the 7 = halo potential. After normalization, 

Ull ~-^ Uw ~- 2C9-i)3 ~ (^Tp (C8) 

with Uiq = —1/4 for (3 = 1. As well known, the scalar virial theorem for a multi-component 
system reads = (|J7n| + //|Wio|)/2, where the interaction energy is W\q = 4ir J r 3 pi(d^o/dr)dr, 
and after normalization, 

/3 2 + 10/? + l 6/?(/3 + l)ln/3 

^ = (^_1)4 + (^_1)5 ' ( C9 ) 

with l^io = -1/10 for (3=1. 



C.2. The 7 = Component 

In this section the normalization constants are the physical scales of the 7 = component, 
i.e., s = r/ro, p = M\/Mq and (3 = n/ro- Consistently with this choice, now = GMq/tq, 
C/n = Mo^n, ^0 = ^00 + V>Iqi , and Jlo = Aqq + pAo\. After normalization and integration one finds 

t 1 + 6s a 20s 3 + 15s 2 + 6s + 1 

700 " 40vr(l + s)6' Ao ° " 80vr(l + s)* " (C10) 

The interaction with the 7 = 1 halo is described by the two contributions 

z 3 1 + s 12s 3 + 6(/? + 5)s 2 -2(/? 2 -8/3-ll)s + /3 3 -5/? 2 + 13/? + 3 

01 " tt((3 - l) 5 (3 + s 4vr(/3- 1) 4 (1 + s) 3 (/3 + s) ' [ ' 

and 

~ W+l) 1 + s 
j4 ° 1 -2vr(/3-l)5 ln ^ + 

6/3Q3 + l)s 3 + 3/3Q3 + 5)Q3 + l)s 2 + (3/3 3 + 25/? 2 + 1(3 + l)s + /3(/3 2 + 10/3 + 1) 

4vr(/3- 1) 4 (1 + s) 3 (/3 + s) ' 1 j 



-25- 



When 13 = 1 

/01 -20vr(l + S )5' ^ " 40vr(l + ' ^ 

The gravitational energy of the 7 = component is given by Uq = Uqo + ^Uqi, where 
Uqo = — 2-7T / pQ^Qr 2 dr is the contribution due to the self-interaction, and Uq\ = —4ir J pQ^\r 2 dr, 
is due to the 7=1 halo potential. After normalization, 

ff - l - TJ - 2/? 2 + 5/3-1 3/3 2 ln/? 

1700 - _ io' Uoi ~~ 203-1)3 + c^iy (C14) 

with C/01 = —1/4 for /3 = 1, obviously equal to f/io- In this case the interaction energy is 
W01 = 47r / r 3 po(d^>i/dr)dr and the virial theorem is = (|C7oo| +/ Lt |W / oi|)/2; after normalization 
one obtains 

17/3 2 + 8/? - 1 3/3 2 (/3 + 3)ln/3 
^ 01 = 209- 1)* 09-l)« ■ ( C15 ) 

For (3 = 1, W 01 = -3/20. 
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Table 1. Critical values of s a = r a /r c for 7 models. The last column refers to the case of a 

dominant central BH. 



7 


NC 


True a 


ssc 


wsc 


WSCbh 





0.354 


0.445 


0.501 


0.556 




1 


0.128 


0.202 


0.250 


0.310 


0.707 


2 


0.000 


0.000 


0.047 


0.107 


0.309 


3 


0.000 


0.000 


0.000 


0.000 


0.000 



Obtained from the model DF. 
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Fig. 1. — The limits on the anisotropy radius for the consistency of the one-component 7 models as 
a function of 7. The solid line represents the lower limit as determined by the NC, and the dotted 
line the upper limit as given by the WSC. For 7 = (0,1,2,3) the black dots connected by the 
dashed line represent the upper limits obtained using the SSC, while the true limits derived from 
the DFs are shown as black squares. The long-dashed line represents the limit on the anisotropy 
radius in presence of a dominant central BH. 

Fig. 2. — The dimensionless DF /i//n for the one-component 7 = 1 model [solid line, Q(0) = 'I'i(O) 
on the abscissae], and in the case of a dominant 7 = halo [lim^oo ^ 3 / 2 /i//n, Q(0) = ^o(O) on 
the abscissae], when (3 = 0.2 (dotted line) and (3 = 3 (dashed line). In the upper panel the globally 
isotropic case is shown, while in the lower panel the anisotropy radius is fixed at r a = O.2677. The 
normalization constant is /n = G~ 3 / 2 M 1 _1//2 7-^' 3 ' /2 . 

Fig. 3. — The dimensionless DF /o//n for the one-component 7 = model [solid line, Q(0) = 
*o(0) on the abscissae, /n = G~ 3 / 2 M 1 ^ 2 r 3 ^ 2 ], and in the case of a dominant 7 = 1 halo 
[lim^ /W _1 /o//N, Q(0) = *i(0) on the abscissae, / N = G- 3 / 2 M[ 1/2 r± 3/2 ], when (3 = 0.2 (dotted 
line) and (3 = b (dashed line). In the upper panel the globally isotropic case is shown, while in the 
lower panel the anisotropy radius is fixed at r a = 0.65ro. 

Fig. 4. — The central value of the dimensionless isotropic [solid line, equation (36)] and anisotropic 
[dotted line, equation (37)] parts of the DF for the halo-dominated 7 = model, as a function of 

(3 = ro/Vi, for (3 > 1. 

Fig. 5. — The minimum value for the anisotropy radius of both components of (1,0) models in the 
case of a dominant halo, as a function of [3 = ro/r\. The anisotropy radius is normalized to the 
scale-length of the specific component. In panel (a) the solid line refers to the 7 = 1 component, 
and the dotted line to the 7 = component. The two dashed lines represent the same quantities, 
when the halo mass is ten times the mass of the specific component. The black dots are the critical 
values of s a for the one-component models. Panel (b) is the enlargement of the small window at the 
top-right of panel (a). The dashed line shows the upper limit on s a due to the halo concentration 
effect. 
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